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SIMULATION OF RADIATION FROM 


LIGHTNING RETURN STROKES: 

THE EFFECTS OF TORTUOSITY 

ABSTRACT 

A Monte Carlo simulation has been developed for the electromagnetio fields 
radiated from a tortuous lightning channel. Tliis was done using a piecewise 
linear model for the chaimel and employing for each element the field radiated 
by a traveling wave on an arbitrarily oriented filament over a conducting plane. 
The simulation reproduces experimental data reasonably well and has been used 
to study the effects of tortuousity on the fields radiated by return strokes. 

Tortuosity can significantly modify the radiated waveform, tending to render 
it less representative of the current pulse and more nearly unipolar than one 
would expect based on the theory for a long straight channel. In the frequency 
domain the effect of tortuosity is an increase in high frequency energy as com- 
pared with an equivalent straight channel. The extent of this increase depends 
on the mean length of the elements comprising the channel and can be significant. 
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SIMULATION OF RADIATION FROM 


LIGHTNING RETURN STROKES; 

THE EFFECTS OF TORTUOSITY 

Introduction 

A Monte Carlo simulation has been developed for the electromagnetic fields 
radiated from a tortuous lightning channel by using a piecewise linear model for 
the chamiel and assuming a traveling wave ("transmission line" model) for the 
current in each element. This simulation has been developed to study the effects 
of tortuosity on the electromagnetic radiation emitted from return strokes in 
cloud-to-ground lightning, and is also proving useful in evaluating the effect of 
other lightning parameters, such as current pulse shape and propagation velocity, 
on the radiated waveforms. 

The simulation is performed by numerically creating a lightning channel out 
of linear segments whose length and orientation are chosen in accordance with a 
specified statistical model. It is assumed that a current pulse with shape Ijj(t) 
propagates along each segment at velocity, v; that is I(t) = I^(t - • r/v). Typi- 

cally this pulse begins at the base of the channel and propagates up toward the 
cloud, possibly dividing and changing velocity at branches. The solution for the 
radiated fields is obtained by summing the fields radiated from each of the linear 
elements, keeping proper track of the phase of the current driving each element 
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as the current pulse moves along the channel. It is assumed that the ground 
plane above which the channel is situated is perfectly conducting (o- = “) so that 
for each element there is also an image element below the plane. The basic 
geometry is illustrated in Figure 1. 

The analytic problem central to the simulation is to determine the electro- 
magnetic fields radiated by the arbitrarilj'^ oriented linear current filaments 
from which the chatmel is constructed. The solution is required in the case of 
a traveling (propagating) current pulse, matiing the problem analogous to cal- 
culating radiation from an arbitrarily oriented traveling wave antenna positioned 
over a conducting plane. Two forms of this solution have been obtained. 


The first solution, and the one used to produce most of the results which 
will be presented here, employs the fraunliofer (i.e. far field) approximation to 
obtain closed form integrals for the fields. This solution is obtained for arbi- 
trary I^(t) and a^rbitrary filament orientation subject only to the constraint that 
the current propagates along the channel. The velocity of propagation as well 
as location of the observer are also arbitrai’y. Tiie time domain form of this 
solution in the special case of an obseiwer located on the conducting surface is 
given below in Equations 1. The general form for the solution and details of its 
derivation are presented in LeVine and Meiieghini (1976a) . 
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where ~ ~ - (L/2c) {'t • Vp^); where p^ is the 

distance from filament center to the observer and L is the length of the fila- 
ment; where the unit vector 't points in the direction of propagation along the 
filament and r^ is the position vector of the center of the filament; and where 
7] = c/v is an index of refraction. The first term in is just the time required 
for the current waveform to propagate halt the length of the filament and the 
second term is the binomial correction to the distance from observer to filament 
center needed to account for the slight differences in location of the end points of 
the filament relative to its center. (See insert in Figure 1.) 


The second solution obtained for the electromagnetic problem has the imique 
attribute of being exact; That is, no mathematical approximateions are necessary 
in arriving at closed form solutions to Maxwell's equations. However, the solu- 
tion applies only to the case where the velocity of propagation, v, of the current 
pulse is equal to c, the velocity of propagation of light in the ambient medium. 
This solution is an extension, to include arbitrary filament orientation and the 
effects of the ground plane, of a solution presented in the time domain by 
Schelkvuioff {Schellcimoff, 1952) and earlier in its frequency domain form by 
Pistollcors (1929) and Beckmann (1931). The form this solution talces in the time 
domain for an observer located on the surface is; 
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and where p , are the distances from the ends of the filament to the obseiwer 
(Figure 1) and r„ ^ are the position vectors which locate the ends of the filament 
mth respect to the reference coordinate system. The general form of the solu- 
tion ajid details of its derivation are to be found in LeVine and Meneghini (1976 a,b) 
This exact solution lias been included as part of the model both to keep track of 
the effect of the approximations made in obtaining the fraunhofer form, and also 
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to give some insight into the near field behaviour of the radiated fields, the re- 
gion where the fraunhofer approximation is inappropriate. 

The current waveform assumed in this simulation is a compound exponential 
model based on the form proposed by Uman (1969) , but modified slightly to be 
continuous at t = 0. The model used here is: 

I = I [e”“‘ - e”^‘] + I, - e“*‘] (4) 

where, for a typical lightning return stroke the parameters are [Uman, 1969]: 

a = 2.0 X 10‘> 

/J = 2.0 X 10® 

V = 1.0 X 10 3 
S = 2.0 X lO'’ 

Iq = 30 ka 
= 2.5 ka 

The first two exponentials in Equation 4 represent the main current pulse in a 
form proposed by Bruce and Golde (1941) with parameters suggested by Deimis 
and Pierce (1964). The third term, Ijexp(-yt), I’epresents intermediate current 
(Uman, 1969); and the remaining exponential, I ^ exp (-St), has been added to 
achieve continuity at t = 0 and doesn't, otherwse, significantly alter the current 
pulse. (See Figure 2 for the shape of the current pulse.) The simulation operates 
independently of tlie actual current waveform, any form being allowed as long as 
it propagates along the channel; however, most of the data obtained to date has 
been based on the preceeding waveform since it appears to reasonably represent 
a composite of reported measurements (Uman, 1969; Dennis and Pierce, 1964). 
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An example of t 5 ^ical data produced by the simulation is shown in Figures 
3 and 4. Figure 3 shows the diree dimensional channel generated by the com- 
puter including, in this case, branches. To generate the channel, the niunber and 
statistics of the main channel elements and branches are specified and then the 
machine produces tlie particular realization by choosing numbers for the length 
and orientation of the elements in accordance with the assigned statistics. The 
number of branches is assigned, but their location is chosen stochastically by 
the computer. The projection of the channel on orthogonal planes is shown in 
the two figures. 

Once the channel has been generated, the computer calculates both the tem- 
poral history of the electromagnetic fields produced by a current pulse propa- 
gating along this chaimel and also the Fourier transform (spectrum) of the radi- 
ated waveform. This is done using the solutions given in Equations I or 2 and 
their frequency domain equivalent, and employing an assigned current waveform 
such as given by Equation 4. The major computational problem is keeping track 
of the proper phase of the current in each of the many elements which may be 
radiating at a given instant. 

An example of the temporal history and spectrum of the electric field wave- 
form radiated from the chaimel shown in Figure 3 and are given in Figure 4, The 
current waveform is as given in Equation 4 with v = c, and the obseiwer is located 
on the y-a.xis 100 km from the origin. Just as the channel appears different to 
observers located at different points so the radiated waveform, measured by 
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observers at different locations is different. However, the shape illustrated in 
Figure 4 is typical of the waveform seen by an observer located on the surface. 
The bottom illustration in Figure 4 represents the spectriun (magnitude of the 
Fourier transform) of the radiated electric field shown at ^he top. The shape 
of the spectrum is typical of an essentially vertical chaimel with the observer 
on the surface; However, the spectrum depends stx-ongly on such parameters as 
channel length and orientation, velocity of propagation, and shape of the current 
pulse. (See LeVine and Meneghini, 1967a for examples and also Figure 9). The 
high frequency asymptote in this example is representative of data obtained with 
a current waveform as described in Equation 4. But at high frequency this data 
is under sampled (because of size limitations in the graph) so that the oscillating 
character evident at high frequency is not quantitatively correct. The number 
and depth of the nulls are not faithfully represented, although the asymptote of the 
pealc value is correctly indicated. The effect of aliasing is especially evident in 
the spectra sho\vn in Figure 5. On the other hand, the first few "oscillations" 
(near 100 Id-Iz) are properly represented. 

The radiated waveform and spectrum presented in Figure 4 have been pro- 
duced using the fraunhofer solution. The exact solution produces results which 
are ostensibly the same in the time domain because the observer is in the far 
field. However, in the frequency domain there is an essential difference at very 
low frequencies. This is because the fraunliofer solution fails for an observer 
closer than several wavelengths from the filament. The fraunliofer solution 
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predicts a decrease in spectrum proportional to frequency in the low frequency 
limit wiiereas the proper behaviour begins to increase as (frequency)"^ at low 
frequency. Thus in Figure 4, the proper spectriun would show an increase at 
frequencies less than about 1 Idlz. For freq\;ency above about 1 Idiz both exact 
and fraunliofer solution agree very well. This difference is illustrated in Fig- 
ure 5 where the spectriun for both the exact and fraunliofer solutions are plotted 
for an observer on the surface and 100 km away from a vertical straight element. 
The current waveform is as given in Equation 4 and the velocity of propagation 
is V - c (necessary for comparison). 
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Comparison with Measured Data: 


The electric field wavefoi'm predicted by the simulation will depend on both 
the particular realization of the channel which happens to be generated, as well 
as on the various parameters with which the simiHation is initialized (i.e. the 
velocity of propagation, current waveform, scale of tortuosity, number and struc- 
ture of branches and location of the observer). In fact, one of the potential uses 
of the simulation is in studying how these various parameters affect the radiated 
waveform. By so doing one, perhaps, can learn to extract channel physics from 
the radiated waveform. On the other hand, a comparison with data can also be 
used as a check on the model, because one would expect that with reasonable 
choices of initial parameters the simulation should produce electric field wave- 
forms which are representative of measured data. 

A comparison of selected simulated and measured electric field waveforms 
is slaown in Figures 6 and 7, The measured data (lower figures) are electric 
field waveforms radiated from first return strokes during cloud-to-ground light- 
ning flashes observed at Keiuiedy Space Center, Florida. These data were ob- 
tained during the Tliunderstorm Research Project (Pierce, 1976) at KSC during 
July 1976 using an electric field measuring system (fast field change meter) 
provided by the University of Arizona (Krider, et al., 1976). The data were 
stored digitally and then reconstructed on im oscilloscope so that photographs 
could be made. The data presented here are representative of first retuim 
strokes observed at KSC. Examples of the simulation are shown at the top of 
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each Figure. The simulated waveform in Figure 6 was produced from the channel 
shown in Figure 3 with the observer at 100 km along the y-axis. The simulated 
waveform shown in Figure 7 was produced from a channel similar to the one 
shown on the bottom of Figure 8. There were no branches in this chamiel and the 
observer was 100 km away from the channel on the X-aXis. However, in this lat- 
ter case, the current waveform was made more peaked by letting a = 6 x 10'* and 
= 10® ; also the velocity of propagation of the current pidse was made to change 
at about 1.2 km such that the velocity was v = c for the first 1. 2 km and v = c/2 
thereafter. The increased structure in the radiated waveform in Figure 7 as 
compared to that in Figure 6 is primarily the result of the modified current 
pulse shape. 

The comparison in Figures 6 and 7 are two representative cases. The com- 
parison is not meant to be statistically significant, which would require large 
samples of unbiased data, but only to show that the simulation can produce qual- 
itatively correct data. Figures 6 and 7 show that assuming reasonable initial 
paranieLers the general features of the simulated waveform are in keeping with 
obsei’ved data. Obviously, many other checks have been made for consistency 
and correctness of the simulation. 
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Effects of Tortuous ity 

The effect of tortuosity on the radiated waveform is especially appropriate 
for study with the simulation because all parameters can be fixed and only the 
channel geometry allowed to vary . This has been done here by creating a sta- 
tistical sample of channels under fixed initial parameters and examining the 
radiation from each. The results are discussed, below first as they appear in 
the time domain and then in the frequency domain. 

A. Time Domain 

The effect of increasing tortuousity on the temporal shape of the electric 
field radiated from a lightning channel is illusirated in Figure 8. The figures 
on the left represent the electric field radiated from the channels shown on the 
right. In each case the current waveform is as given in Equation 4 and the pulse 
propagates along the channel at velocity v = c/2. The observer is located on the 
X-axis 100 km from the channel. 

The top row of figures are the results due to a long straight vertical channel 
(i.e. no tortuosity). For such a channel the radiated fields mathematically sepa- 
rate into two terms (see Equations 1), one of which corresponds to radiation 
emanating from the base of the channel and beginning when the pulse starts to 
propagate up the channel, and the second corresponding to radiation from the 
top of the channel but not beginning until the pulse has reached the top. The 
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wave form radiated from each end point is proportional to the shape of the cur- 
rent pulse. Thus, initially, before the pulse reaches the top of the channel, the 
shape of the I'adiated waveform is the same as the current pulse because only 
one term is contributing. The time (approximately 40 ^ s) at which tlie current 
pulse reaches the top of the chaimel is quite evident in the upper figure. At this 
time the second radiation pulse begins (with identical shape but opposite sign); 
and the subsequent electric field waveform is their sum. 

Wlien the chaimel is tortuous (lower two figiu'es), each linear element be- 
haves as described above. The total field at the observer is then the sum of all 
elements radiating at a given time. Matliematically this is equivalent to many 
"point" radiators located at the channel icinks (i.e. points at which the linear 
segments join). Wlien the adjoining segments are nearly parallel, the radiation 
from the top of one nearly cancels that from the bottom of the other, and tlie net 
effect is small. This is seen in the "low" tortuousity chaimel shown on the middle 
row of graphs in Figure 8. However, when the angle between adjoining segments 
is lai'ge as can happen in highly tortuous channels, the radiation from eacli jimc- 
tion can be significant. The net effect, shown for the higlily toxtuous channel at 
the bottom of Figui'e S, can alter the basic shape of the radiated waveform as 
inferred from the long straight channel. The radiated waveform loses its simi- 
larity to the current pulse and the pronomiced "turaing off" at the top of the 
channel which is apparent if the top figure becomes less significant. The proc- 
ess of extracting information about the shape of the cimi'ent pulse from the 
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the radiated waveform will become more and more difficult as tortuousity in- 
creases. Even the leading edge can be altered, rendering difficult the associa- 
tion of rise time as measured from the radiated waveform with that of the cur- 
rent pulse. They are the same imtil the first "kink" is reached, but not so 
thereafter. 

B. Frequency Domain 

The effects of tortuosity as discussed above in terms of the temporal shape 
of the radiated waveform, are also manifested in the frequency domain. Two 
factors dominate the effect in the frequency domain: filament orientation and 
filament length. The tortuous channel, as compared to the vertical straight 
channel, acquires an increasingly greater proportion of horizontal orientation 
as the tortuousity is increased due to the random orientation of the component 
elements. Also the more variability per unit length in a given chamiel, the 
shorter the elements must be to accommodate the frequent changes in direction. 

Orientation of the channel can be an important factor determining the shape 
of the spectrum of the radiated electromagnetic fields. This is illustrated in 
Figure 9 where the spectrum of the electric field radiated from a straight (i.e. 
single element) vertical channel and straight horizontal channel are shown. In 
each case the current waveform is as given in Equations 4 with v = c, a channel 
length of 6 km (the horizontal channel is 4 km above the surface) and the ob- 
server 100 Ion away. The horizontal channel has relatively more energy at high 
frequency and less at low frequencies as compared with an otherwise identical 
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vertical channel. Since the tortuous channel can be thought of as adding hori- 
zontal components to what otherwise would be a long strtiight vertical channel, 
one would expect to see a shift of energy out of the lower frequencies and toward 
higher frequency in the tortuous channel. It should be noted ho^vever, that this 
oi’ientational effect is velocity dependent and is much less pronounced at veloc- 
ities less than about v = c/2 than shown in Figiu’e 9. 

An increase in high frequency energy is also caused by shortening the length 
of a given element. The spectrum scales roughly as frequency multiplied by ele- 
ment length {LeVine and Meneghini, 1976a) so that for a given orientation the 
spectrum shifts toward high frequency as the element length is decreased. In a 
tortuous chamiel one would expect both smaller mean element length and an in- 
crease in relative horizontal orientation. The net effect due to both is a shift of 
energy toward higher frequencies. At velocities of propagation much less than 
V = c, element length will be the dominant effect. 

This increase in high frequency energj'^ is illustrated in Figure 10 where the 
spectrum of the electric field radiated from a vertical straight channel is com- 
pared with that from an otherwise identical tortuous channel. At frequencies 
greater than a few hiuidred kHz the spectrum associated with the tortuous chan- 
nel lies above that due to the straight vertical channel. This shift of energy 
toward higli frequencies can be significant, an order of magnitude or two easily 
being produced with reasonable lightning parameters. 
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01 course, what has just been described in the frequency domain is only 
aiiothoi* way of looking at the behaviour of the radiated fields in tlie time domain. 
Thus, in the time domain the effect of tortuousity is an inci'eased fluctuation of 
the radiated signal (per unit time) . This is due to the variability in the channel 
ovlhutatign aiid decrease in the length of the elements. Signals with a high de- 
gree of variability require more high fi'equencj' content for their reproduction 
than do more smoothly varying signals. 
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Summary 

A Monte Carlo simulation has been developed to model radiation from light- 
ning return strokes. The simulation employs a piecewise linear model for the 
chaimel and assumes a current pulse which propagates along the channel. The 
results of the simulation compare favorably with measured data. 

The simulation has been used to study the effects of tortuosity on the elec- 
tric field waveforms radiated from return strokes. Tortuousity tends to make 
the radiated fields less representative of the current pulse and more nearly uni- 
polar than one would expect from calculations using a straight channel (i.e. single 
element). The effect of tortuousity in the frequency domain is an increase in 
high frequency energy and reflects the increased structure (fluctuations) in the 
time domain. The magnitude of this increase depends on the mean length of the 
elements comprising the chamiel. 
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Figure 4. An example of the electric field waveform 
produced by the simulation in the time domain and 
frequency domain. The channel used is that shown in 
Figure 3 nnd the current waveform is as given in Fig- 
ure 3. The velocity of propagation of the current 
pulse is V = c and the observer is on the y-axis, 100 
km away. 
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Fijiurc <>. Coiiipiu'ison of simulatetl electric field waveform (top) 
with a measured waveform (bottom). The measured data 
rei)resents radiatiou from a first return stroke at KSC. Florida. 
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Figure 7. Comparison of simulated electric field waveform (top) 
with measured waveform (bottom). The measured data are 
from first return strokes observed at KSC, Florida, 
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Figure 8, The effects of increasing tortuousity on the 
radiated electric field waveform. 
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Figure 9. The effects of channel orientation on the spectrum of the radiated electric field. 
The observer is on the surface 100 km from the origin in each case. Each channel is 6 km 
long and the horizontal channel is 4 km above the surface. 
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Figure 10. The effects of tortuousity in the frequency domain. Thi. electric field radiated 
from the tortuous channel has more high frequency content than a similar straight channel. 



